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We observe evidence for D°-D° mixing by measuring the difference in apparent lifetime when a 
D° meson decays to the CP eigenstates K + K~ and tv + ty~, and when it decays to the final state 
K~tt + . We find y C p = (1.31 ± 0.32(stat.) ± 0.25(syst.))%, 3.2 standard deviations from zero. We 
also search for a CP asymmetry between D° and D° decays; no evidence for CP violation is found. 
These results are based on 540 fb _1 of data recorded by the Belle detector at the KEKB e + e~ 
collider. 

PACS numbers: 13.25.Ft, 11.30.Er, 12.15.Ff 



The phenomenon of mixing between a particle and its 
anti-particle has been observed in several systems of neu- 
tral mesons 0, Q: neutral kaons, B^, and most recently 
B® mesons. In this paper we present evidence for D°-D° 
mixing [3j. 

The time evolution of a D or D° is governed by the 
mixing parameters x — (Mi — Mzj/Y and y = {T\ — 
T2)/2T, where M\$ an d are the masses and widths, 
respectively, of the mass eigenstates, and V = (ri+r 2 )/2. 
For no mixing, x = y = 0. Within the Standard Model 
(SM) the rate of D-mixing is expected to be small due 
to the near degeneracy of the s and d quark masses rel- 
ative to the W mass, and the small value of the b quark 
couplings. Predictions for x and y are dominated by 
non-perturbative processes that are difficult to calculate 
dH- The largest predictions are |at|, \y\ ~ 0(1O~ 2 ) 
Loop diagrams including new, as-yet-unobserved parti- 
cles could significantly affect the experimental values Q . 
CP-violating effects in D-mixing would be a clear signal 
of new physics, as CP violation (CPV) is expected to be 
very small in the SM 0]. 

Both scmilcptonic and hadronic D decays have been 
used to constrain x and y [l| . Here we study the decays to 
CP eigenstates D° — > K + K~ and D° — > 7r + 7r~; treating 
the decay time distributions as exponential, we measure 



the quantity 



t(K-tt+) 
VCP ~ t{K+K-) 1; 



(1) 



where t(K + K ) and t(K ir + ) are the lifetimes of D° —> 
K+K- (or TT+n-) and D° -> K-jt + decays 0. It can 
be shown that yep — y cos <f> — ^Amx sin tf> [9( , where Am 
parameterizes CPV in mixing and is a weak phase. If 
CP is conserved, Am = = and yep = y- To date 
several measurements of yep have been reported [l(J; the 
average value is ^2 standard deviations (cr) above zero. 
Our measurement yields a nonzero value of yep with 
> 3a significance. We also search for CPV by measuring 
the quantity 



Ar = 



t(D° -> K~K+) - t{D° -> K+K- 



-(D° 



r(D° 



K+K- 



(2) 



this observable equals Ar = ^AmV cos <p — ismi 

Our results are based on 540 fb _1 of data recorded 



by the Belle experiment at the KEKB asymmetric- 
energy e + e~ collider [l2l ]. running at the center-of-mass 
(CM) energy of the T(45 l ) resonance and 60 MeV below. 
To avoid bias, details of the analysis procedure were final- 
ized without consulting quantities sensitive to yep and 
A T . 
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FIG. 1: M distribution of selected events (with |Ag| < 
0.80 MeV and a t < 370 fs) for (a) K + K~, (b) K~ir + and 
(c) tt + 7v~ final states. The histogram shows the tuned MC 
distribution, (d) q distribution (with \AM\/ani < 2.3 and 
a t < 370 fs) for the K + K~ final state, (e) Normalized dis- 
tribution of errors a t on the decay time t for D° — > K~ir + , 
showing the construction of the resolution function using the 
fraction ft in the bin with at = (f) Fitted lifetime of 
D° mesons in the K~ir + final state in four running periods 
with slightly different conditions, and the result of a fit to a 
constant. The world average value (W.A.) is also shown. 



The Belle detector is described in detail elsewhere 11 1 : 
it includes in particular a silicon vertex detector 13j |. 
a central drift chamber, an array of aerogel Chercnkov 
counters, and time-of-flight scintillation counters. We 
reconstruct D* + — » D°ir+ decays with a characteristic 
slow pion 7r s , and D° — > K + K~, K~~tt + , and 7r + 7r~. 
Each track is required to have at least two associated 
vertex detector hits in each of the two measuring coor- 
dinates. To select pion and kaon candidates, we impose 
standard particle identification criteria 14|. D° daugh- 
ter tracks are refitted to a common vertex, and the D° 
production vertex is found by constraining its momen- 
tum vector and the ir s track to originate from the e + e~ 
interaction region; confidence levels exceeding 10 -3 are 
required for both fits. A D* momentum greater than 
2.5GeV/c (in the CM) is required to reject D-mesons 
produced in £?-meson decays and to suppress combinato- 
rial background. The proper decay time of the D° can- 
didate is then calculated from the projection of the vec- 
tor joining the two vertices, L, onto the D° momentum 
vector, t — mjjaL ■ p/p 2 , where mjja is the nominal D° 



mass. The decay time uncertainty at is evaluated event- 
by-event from the covariance matrices of the production 
and decay vertices. 

Candidate D° mesons are selected using two kinematic 
observables: the invariant mass of the D° decay prod- 
ucts, M, and the energy released in the D* + decay, 
q = {Md* — M — m^)c 2 . Mjj* is the invariant mass 
of the D°ir s combination and m n is the ir + mass. 

According to Monte Carlo (MC) simulated distribu- 
tions of t, M, and q, background events fall into four 
categories: (1) combinatorial, with zero apparent life- 
time; (2) true D° mesons combined with random slow 
pions (this has the same apparent lifetime as the signal); 
(3) D° decays to three or more particles, and (4) other 
charm hadron decays. The apparent lifetime of the latter 
two categories is 10-30% larger than r^o. Since we find 
differences in M and q distributions between MC and 
data events, we perform fits to data distributions to ob- 
tain scaling factors for the background fractions and sig- 
nal widths, and then tune the signal fractions and shapes 
in the MC event-by-event. 

The sample of events for the lifetime measurements 
is selected using \AM\/<jm, where AM = M — muo; 
\Aq\ = q — (m£)*+ — m D a — m v )c 2 ; and at- The invariant 
mass resolution <tm varies from 5.5-6.8 MeV/c 2 , depend- 
ing on the decay channel. Selection criteria are chosen to 
minimize the expected statistical error on yep, using the 
tuned MC: we require \AM\/a M < 2.3, \Aq\ < 0.80 MeV, 
and at < 370 fs. The data distributions and agreement 
with the tuned MC are shown in FigsQJa)-(d). We find 
111 x 10 3 K+K-, 1.22 x 10 6 K~ir+, and 49 x 10 3 tt+tt" 
signal events, with purities of 98%, 99%, and 92% respec- 
tively. 

The relative lifetime difference yep is determined from 
D° — » K + K~, K~ir + , and ir + ir~ decay time distri- 
butions by performing a simultaneous binned maximum 
likelihood fit to the three samples. Each distribution is 
assumed to be a sum of signal and background contribu- 
tions, with the signal contribution being a convolution of 
an exponential and a detector resolution function, 



dN/dt = ^is 



r*'/ T • R(t - 1') dt' + B{t), (3) 



The resolution function R(t — t') is constructed from the 
normalized distribution of the decay time uncertainties 
a t (see Fig. HJe)). The a t of a reconstructed event ide- 
ally represents an uncertainty with a Gaussian probabil- 
ity density: in this case, we take bin i in the at distri- 
bution to correspond to a Gaussian resolution term of 
width crj , with a weight given by the fraction fi of events 
in that bin. However, the distribution of "pulls", i.e. 
the normalized residuals (t rec — t gcn )/at (where i rcc and 
tgen are reconstructed and generated MC decay times), 
is not well-described by a Gaussian. We find that this 
distribution can be fitted with a sum of three Gaussians 
of different widths cF uU and fractions Wk , constrained to 
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FIG. 2: Results of the simultaneous fit to decay time dis- 
tributions of (a) D° -> K+K~, (b) D° -> K~ir+ and (c) 



D 



decays. The cross-hatched area represents back- 



ground contributions, the shape of which was fitted using M 
sideband events, (d) Ratio of decay time distributions be- 
tween D° -> K + K~,tt + -k' and D° -> K~n + decays. The 
solid line is a fit to the data points. 



the same mean. We therefore choose a parameterization 



R(t-t') = ^2fi^2w k G{t-t';a ik ,to), 



(4) 



fe=i 



with <7ik = Sfeffe fj, where the s k are three scale factors 
introduced to account for differences between the simu- 
lated and real c£ u11 , and to allows for a (common) offset 
of the Gaussian terms from zero. 

The background B{t) is parameterized assuming two 
lifetime components: an exponential and a 8 function, 
each convolved with corresponding resolution functions 
as parameterized by Eq. ((4]). Separate B(t) parameters 
for each final state are determined by fits to the t dis- 
tributions of events in M sidebands. The tuned MC is 
used to select the sideband region that best reproduces 
the timing distribution of background events in the sig- 
nal region. We find good agreement between the tuned 
MC and data sidebands, with a normalized x 2 of 0.85, 
0.83 and 0.83 for KK, Kir, and itit respectively. 

The R(t — t') and background parameterizations are 
validated using MC and the large D° — ► K~ir + sam- 
ple selected from data. In the simulation, the ratio of 
scale factors s k (k = 1,2,3) is consistent between decay 
modes, within small statistical uncertainties. The offset 
to is also independent of the final state, but it changes 
slightly for simulated samples describing different run- 
ning periods. Four such periods, coinciding with changes 



to the detector, have been identified based on small vari- 
ations of the mean t value for D° — * K~ir + in the data. 
We perform a separate fit to each period and average the 
results to obtain the final value of yep ■ The free parame- 
ters of each simultaneous fit are: r^o, yep, three factors 
Sk for the K~ir + mode and two terms that rescale the 
s k in the K + K~ and tt + it~ channels, the offset to, and 
normalization terms for the three decay modes. Fits to 
the D° — > K~n + sample show good agreement with the 
parameters of R(t — t') obtained from simulation. 

For the second running period, we modify Eq. ^ to 
add mode-dependent offsets At between the first two 
Gaussian terms, making the resolution function asym- 
metric; these three parameters are also left free in the 
fit. We find that such a function is required to yield 
the D° — ► K~ir + lifetime consistent with that in the 
other running periods. (This behaviour has been repro- 
duced with a MC model including a small relative mis- 
alignment of the vertex detector and the drift chamber.) 
The lifetime fit results are shown in Fig. [TJf ) : the mean, 
t d o = (408.7 ± 0.6 (stat.)) fs, is in good agreement with 
the current world average, (410.1 ± 1.5) fs Q. 

Fits to the D° -> K+K~, A'~7r+ and tt+tt" data 
for the four running periods are shown in Fig. [2]Ja)- 
(c), by summing both the data points and the fit func- 
tions. Averaging the fit results, we find yep — (1.31 ± 
0.32 (stat.))%, 4.1 standard deviations from zero. The 
agreement between the data and the fit functions is 
good: x 2 / n dof = 1-08 for ridof = 289 degrees of free- 
dom. Fitting K + K~ / K~tt + and iy + tt~ / K~tt + events 
separately we obtain yep = (1.25 ± 0.39 (stat.))% and 
yep = (1.44 ± 0.57 (stat.))% respectively, in agreement 
with each other. The yep values for the four running 
periods are also consistent, with x 2 / n dof — 1-53/3. 

To measure the CPV parameter Ay we separately de- 
termine the apparent lifetimes of D° and D° in decays 
to the CP eigenstates; the data is fit in four running 
periods as for yep- To ensure convergence of the fits, 
despite the much smaller event sample, the scale fac- 
tor for the widest Gaussian S3 is fixed to the value ob- 
tained from the yep fit in each case. We obtain Ar — 
(0.01 ± 0.30 (stat.))%, consistent with zero; the quality 
of the fit is good, with x 2 / n dof = 1-00 for ridof = 390. 
Separate fits to the two CP eigenstates find compati- 
ble values: A T = (0.15 ± 0.35 (stat.))% for K+R- and 
-(0.28 ± 0.52 (stat.))% for vr+Tr". 

The behaviour of the fits has been tested in various 
ways using MC simulation. Fits to signal events simu- 
lated with yep = reproduce this value (and the gener- 
ated Tdo) even for a sample much larger than the data, 
with (x 2 1 ndof , ndof) = (1.11,285). Using samples of the 
same size as the data, with background included, we find 
a satisfactory fit, (x 2 / n dof,nd f) = (1.18,289), with a 
statistical uncertainty in agreement with the error from 
the fit to the data. Results obtained on reweighted MC 
samples that cover a wide range of yc p values agree with 
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the input within ±0.04%. 

The effect of the resolution function on the measured 
Ucp has been tested by replacing the parameterization 
in Eq. Q with a single Gaussian. This describes the 
data poorly and leads to a 3.9% shift in the fitted Tjjo 
for a simulated D° — > K~tt + sample; however, the cor- 
responding shift in yep is only 0.01%. This shows that 
the yep value returned by the fit is robust against im- 
perfections in the parameterization of R(t — t'). 

The estimated systematic uncertainties are summa- 
rized in Table |TJ We test for acceptance variations with 
decay time by fitting the generated decay times of recon- 
structed MC events. We find no deviation, but conserva- 
tively assign the MC statistical error on y C p (±0.12%) 
to this source. Another contribution is due to the choice 
of equal to offsets in different decay modes: relaxing this 
assumption leads to yep changes of ±0.14%. Variation of 
the D° mass windows changes yep by less than ±0.04%. 
The effect of differences between backgrounds in the sig- 
nal and sideband regions is studied by repeating the fits 
using MC backgrounds from signal regions; small shifts 
in the data sidebands used to determine B(t) are also 
made. The largest resulting change in yep, ±0.09%, is 
quoted as the systematic error due to the background 
description. Potential correlations between apparent life- 
times and opening angle distributions (which differ be- 
tween modes) have a small effect on yep- ±0.02%. 

The uncertainty due to the finite number of sideband 
events, ±0.07%, is estimated by varying bin contents ac- 
cording to Poisson statistics and repeating the fits. Com- 
paring alternative fits where all running periods use the 
symmetric resolution function Q, and the asymmetric 
function presently used for the second period, we assign 
an additional uncertainty of ±0.01%. Varying selection 
criteria produces observable effects only in high statistics 
MC samples, in the a t and \AM\/cru cases. The result- 
ing ±0.11% changes in yep are conservatively assigned 
as systematic errors. Finally, varying the binning of the 
decay-time distribution produces a small effect, ±0.01%. 
Adding all terms in quadrature, we obtain a systematic 
uncertainty on yep of ±0.25%. The same sources domi- 
nate for Ar, but yield a smaller total systematic uncer- 
tainty, ±0.15%. 

In summary, we measure the relative difference of the 
apparent lifetime of D° mesons between decays to CP- 
even eigenstates and the K~tt + final state to be 

ycp = (1.31 ± 0.32(stat.) ± 0.25(syst.))%. (5) 

Combining the errors in quadrature, we find a confidence 
level of only 6 x 10~ 4 for the yep = hypothesis. Wc 
interpret this result as evidence for mixing in the D°- 
D° system, regardless of possible CPU. The effect is 
presented visually in Fig. [2^d), which shows the ratio 
of decay time distributions for D° — » K + K~ , ir + ir~ and 
D° — > K~n + decays. We also search for CP violation by 
separately measuring decay times of D° and D° mesons 



TABLE I: Sources of the systematic uncertainty for yep and 



Source 


Aycp [%] 


AA r [%] 


acceptance 


0.12 


0.07 


equal to 


0.14 


0.08 


M window position 


0.04 


< 0.01 


signal/sideband background differences 


0.09 


0.06 


opening angle distributions 


0.02 




background distribution B(t) 


0.07 


0.07 


(a) symmetric resolution function 


0.01 


0.01 


selection variation 


0.11 


0.05 


binning of t distribution 


0.01 


0.01 


Total 


0.25 


0.15 



in CP-even final states. We find an asymmetry consistent 
with zero, 

A T = (0.01 ± 0.30(stat.) ± 0.15(syst.))%. (6) 
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